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Abstract 

Analytical expressions for the dispersion relations and polarizations of low-frequency waves in 
magnetized plasmas based on two-fluid model are obtained. The properties of waves propagating 
at different angles (to the ambient magnetic field Bq) and (3 (the ratio of the plasma to magnetic 
pressures) values are investigated. It is shown that two linearly polarized waves, namely the fast 
and Alfven modes in the low-/? (/? <C 1) plasmas, the fast and slow modes in the /? ~ 1 plasmas, 
and the Alfven and slow modes in the high-/? (/? S> 1) plasmas, become circularly polarized at 
the near-parallel (to Bq) propagation. The negative magnetic-helicity of the Alfven mode occurs 
only at small or moderate angles in the low-/? plasmas, and the ion cross-helicity of the slow 
mode is nearly the same as that of the Alfven mode in the high-/? plasmas. It also shown the 
electric polarization 5Ez/5Ey decreases with the temperature ratio Tg/Tj for the long-wavelength 
waves, and the transition between left- and right-hand polarizations of the Alfven mode in Tg/Tj A 
0 plasmas can disappear when Tg/Tj = 0. The approximate dispersion relations in the near¬ 
perpendicular propagation, low-/?, and high-/? limits can quite accurately describe the three modes. 
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I. INTRODUCTION 


It is well known that in homogeneous magnetized two-fluid plasmas three electromag¬ 
netic modes with frequency less than the electron cyclotron frequency exist These 

include the fast, Alfven (or intermediate) and slow modes, according to their different phase 
velocities j^, Q . Their dispersion relations can be obtained from a general one based on the 
Hall-MHD model [I|-[5|. Recently, 0 obtained the same relation from a simpler formula¬ 
tion involving a two-dimensional current density vector. The general dispersion relation for 
even lower frequency modes (with wave frequency u less than the ion cyclotron frequency 
(Uci) have also been derived using different formulations However, a comprehensive 

investigation of the wave polarizations is still lacking. 

In this paper we present analytical expressions of the dispersion relations and polariza¬ 
tions using an approach similar to that of Ref. jsl . We shall consider the dense-plasma limit 
y\l(? -C 1, where Va is the Alfven speed and c is the light speed, so that the displacement 
current in the Ampere’s law can be ignored The resulting analytical expressions are 

useful for analyzing the properties of low-frequency waves in different plasmas. 

In the next section we present the derivation of the dispersion relations and polariza¬ 
tions of the waves. In Sec. Ill the properties of the waves propagating at different angles 
and different (3 regimes are discussed. The main results are summarized in Sec. IV. The 
Appendix gives the approximate dispersion relations in the near-perpendicular propagation, 
I 0 W-/S (/9 -C 1), and high-/? (/? S> 1) limits, where /? is the ratio of the plasma to magnetic 
pressures. 


II. DISPERSION RELATIONS AND POLARIZATIONS 


Linearized two-fluid and Maxwell’s equations are 


manodtS\a = noQa ((5E -F 6\a X Bo) - V6Pa, (1) 

dtSua = -V ■ (rio^Vo), (2) 

V X hB = /iohJ, (3) 

V X 5E = -dt6B, (4) 


where the subscripts a = i,e denote ions and electrons, respectively, ma is the mass, Qq 
is the charge, 6Pa = K'jaTaSna is the thermal pressure, k is Boltzmann constant, is 
the temperature, Sua is the perturbed number density, is the perturbed velocity, 5J is 
the perturbed current density, (5E and 5B are the perturbed electric and magnetic helds, 
respectively, Bq = Rqz is the ambient magnetic held, and no is the ambient number density. 
As mentioned, the displacement current is neglected. The quasi-neutrality condition (fn* = 
drie = 5n shall also be used. In the study an electron-proton plasma is considered, namely 
qe = —e and g* = e. 

We shall consider plane waves, so that 6f oc 6fkexp{—iut -|- ?k ■ r), where oj is the wave 
frequency and k = {k±ex + k^Gz) is the wave vector. We can obtain from Eq. ([1]) the 
perpendicular and parallel (to Bqz) huid velocities 
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The current density (5J = rioe (5vi — 5ve) can then be expressed as 
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where Q = me/rrii, Ao = 1 - A 2 = 1 - + 7 eTe, Ti = '^{Ti/Tt and 

Tg = 'jeTe/Tt. Combining Eqs. (3) and (4) leads to 
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From Eqs. ([7]) - (E]), we get for the electric held and number density perturbation, 
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where Ai = 1 — /oj‘1^ and other dehnitions are: 
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Inserting above electric held components into the number density equation that is derived 
from Eqs. ([2]), ([5]) and (E]), 


[(Ao + plkl) - AoE^Afc^] — = + "^SEy + tk,—Ao5E,, (14) 

no EoUci Bq rrii 

the general dispersion relation can be expressed as 
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where = pj + pi, Pi = VTi/'^d is the ion gyroradius, ps = VtsXci is the ion acoustic 
gyroradius, A, is the ion inertial length, Ag is the electron inertial length, Vr* = \fXXiXi 
is the ion thermal speed, Vts = y/^rfeX/X is the ion acoustic speed, hr = i/nT^/n^is 


the sound speed, and /3 = V^/V|. With respect to the existing ones [i|, [i, 0,0,0, li, 13 


Eq. fllSp represents a more general description of the low-frequency electromagnetic waves. 
Three roots for correspond to the fast {j = 0), Alfven [j = 1), and slow {j = 2) modes, 

or 0[I3, 


ujI = 2p^Xos Qcos 1 - yj) + ^, J = 0,1, 2, (16) 

with p= {B^ - SAC) /(9A2) and q = {9ABC - 2B^ - 27A^D) /{bAA^). If we set fc — )■ 00, 
Eq. flT5]) yields two resonances {ck/u —>■ cx)): the ion cyclotron resonance u = oJdCosO and 
the electron cyclotron resonance u = \u:ce\ cos6. 

If we neglect the electron inertial terms terms of the order of Q, Eq. flT5|l 

recovers the Hall-MHD dispersion relation [l2|, where only the ion cyclotron resonance 
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exists. For the high oblique propagation, low-/9 and high-/9 limits, the approximate dispersion 
relations of the three modes are given in the Appendix. Eq. ([TSjjcan also be reduced to the 
well-known results in the cold two-fluid plasmas (Tj = Tg = 0) [l3|, [14 , 


Once the electric held perturbation 0131) and the dispersion relation (IT^ are known, the 
magnetic held and velocity perturbations can be also expressed in terms of the number 
density perturbation, 
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Note that we can explore the linear relation between arbitrary two variables through the 
eigenfunctions flT^ and ([T7|) — (IT^ . For example, the polarizations of electromagnetic helds 
are 
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Parallel Waves 


At parallel propagation, k = fc^z, Eq. flT5|l is written as 
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and ion aconstic wave 


= V^kH (1 + Q). (24) 

Note that the dispersion relation fl2^ can be directly derived from Eqs. flTOl) and flTTl) : fl24l) 
can be derived by nse of Eqs. flT^ and ffT4D . 

The left- and right-hand waves have the perpendicnlar pertnrbations 
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whereas the ion aconstic wave has the parallel pertnrbations 
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B. Perpendicular waves 


When the wave propagates at the perpendicnlar direction, k = fc_LX, only one mode exists 

2 _ 1 + Q + (1 + Q + Agfc 5_)/3 2k 2 l'07'i 

Its polarization properties are 
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III. DISCUSSION 

At the parallel propagation, the ion aconstic wave can interact with the right/left 
circnlarly-polarized waves at interaction points where their u/ud and Xik, are eqnal as 
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FIG. 1: Wave frequency and phase velocity of the parallel right-hand circularly-polarized wave 
(solid lines), left-hand circularly-polarized wave (dashed lines), and ion acoustic wave (dotted 
lines) in the plasmas with Tj = Tg but different /?: /3 = 10“^, 1 and 10^. 


shown in Fig. (1). A mode transition can occur at the interaction point. The mode transi¬ 
tion can happen among three oblique waves in Figs. (2) —(4). Also, Figs. (2) —(4) include 
the wave electromagnetic polarizations as well as the magnetic helicity a and the ion cross- 
helicity aci 

A; (A ■ hB*) 2cos6 {idBx/6By) 

^ ~ 552 “ cos^e + {5BI/5BI) ’ 

and 

2Re (5vj ■ 5v^) 
dvj + 5v% 


where A denotes the vector potential and Svb = dB j is the magnetic held pertur¬ 
bation in the velocity unit. 

Fig. (2) presents the dispersion relations and polarizations of the three oblique waves at 
different angles in the low-/? plasmas where /? = 10“2 and T* = Tg. It shows that approximate 
dispersion relations (1A4|1 and (1A6|1 can describe the exact one (ITHll well. 

The fast mode corresponds to the fast magnetosonic wave as w <C ojd and the whistler 
wave as Ud ^ OJ ^ IcUggl cos 0. At the electron cyclotron frequency u = IcUggl cos 6*, the 
fast mode is the electron cyclotron wave. Furthermore, the electron cyclotron wave can 
change to a (quasi-) electroacoustic wave extending to higher frequency u > IcUggl cos 6* j^. 
It is interesting to see that the fast magnetosonic wave has {\dEy/dEx \ , \dBx/dBy\) > 1 
as Xik < 10“^, and {\dEy/dEx\ ,\dBx/dBy\) ~ 1 as 10“^ < Xik -C 1 at the near-parallel 
propagation. \dE/dB\ ~ Va and a ~ 0 for the fast magnetosonic wave; \dE/dB\ > Va and 
cr ~ 1 for the whistler and electron cyclotron waves. In addition, aci — — cos 6 for the fast 
magnetosonic wave [^, while aci — 0 for the whistler and (quasi-) electroacoustic waves. 

The Alfven mode is the shear Alfven wave at a; -C Ud and the (quasi-) electroacoustic 
wave at ujd < OJ < |ci;cg|cos6* until a transition into the electron cyclotron wave at a; = 
|a;cg|cos6'. At near-parallel 9 = 5° and oblique 9 = 45° cases, the phase velocity of the 
(quasi-) electroacoustic wave is about the sound speed. It is the Alfven speed at the high 


111 ] called the high oblique mode at Ud < oj < IcUgglcos^ 


oblique angle, therefore. Ref. 
as the shear Alfven wave. Note that an ion cyclotron wave a; ~ cjg 


arises at near-parallel 
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FIG. 2: Wave dispersion relation and polarization at three propagating angles {9 = 5°, 45° and 
85°) in the low-/3 plasmas with /3 = 10“^ and Tg = Tj, where the solid, dashed and dotted lines 
denote the fast, Alfven and slow modes, respectively. Thin solid lines in Panels (al), (bl) and (cl) 
correspond to approximate dispersion relations ()A4I) and ()A6I) in the low-/3 (/3 ^ 1) limit. 


propagation (Panel (al)). Electromagnetic polarizations are \5Ex\ 3> \5Ey\, \^By\ ^ \5Bx\ 
and \5E/5B\ = I /4 at a; <C ojd] at < uj < \ujce\ cos9, \6E/6B\ becomes much larger than 
Va- At 6 = 5° and 45°, the magnetic-helicity a decreases firstly from a = 0 to a ~ —1 at 
oj <C Uci, and then increases with increasing Aj/c at < uj < \oJce \ cos 9] at 0 = 85°, <7 = 0 
is nearly unchanged at a; <C Ud, and it becomes increasing at a; > Ud until reaching cr = 1 
corresponding to the electron cyclotron wave. Besides, the ion cross helicity aci depends on 

8 
















































































0=5 


0=45 


0 = 85 " 



FIG. 3: Wave dispersion relation and polarization at three propagating angles {9 = 5°, 45° and 
85°) in /3 = 1 plasmas with Tg = Tj, where the solid, dashed and dotted lines denote the fast, 
Alfven and slow modes, respectively. Thin solid lines in Panels (al), (bl) and (cl) correspond to 
the approximate dispersion relations (IA2|) and (jASj). 

the wave scale, e.g., aci = —1 as Aj/c 1 and aci = 0 as Ajfc 1. 

The slow mode corresponds to the slow magnetosonic wave at ca Ud, where \5Ex\ < 
\5Ey\, \5Bx\ > \SB„\, \6E/6B\ ~ Vr and <7 = 0. It tnrns to the ion cyclotron wave at 
u = UciCos6 [lil, where \6Ex\ S> \SEy\, \6Bx\ ~ \SBy\, \6E/6B\ > Va and a = —1. At 
6 = 5°, the electric polarization \6Ey/6Ex\ has an increment at the transition position where 
the slow magnetosonic wave changes to the ion cyclotron wave; however, there is no such 
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FIG. 4: Wave dispersion relation and polarization at three propagating angles {9 = 5°, 45° and 
85°) in the high-/3 plasmas with (3 = 10^ and Tg = Ti, where the solid, dashed and dotted lines 
denote the fast, Alfven and slow modes, respectively. Thin solid lines in Panels (al), (bl) and (cl) 
represent the analytical dispersion relations ()A2p and (|A3p corresponding to the high-/3 (/3 ^ 1) 
limit. 

increment a.t 6 = 45° and 6 = 85° cases. 

Figs. (3) and (4) present the dispersion relations and polarizations in /3 = 1 and high-/? 
(/? = 10^) plasmas. Here the Alfven mode interacts with the fast mode only. Although the 
validity condition for approximate dispersion relations flA2p and (|A3|) is kl/k^, » 1 or 
(3^1, these expressions can describe the wave dispersion relations at /? = 1 as shown in 
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Fig. (3). 

Several mode properties in Fig. (3) are obviously different from that in the low-/? plasmas 
(Fig. (2)). For example, to the near-parallel waves at Q.loJd < oj < ojd, two circularly 
polarized {\SEx\ ~ \SEy\) modes are the fast and slow modes in /? = 1 plasmas but the fast 
and Alfven modes in the low-/? plasmas. Here the fast (slow) mode exhibits the right-hand 
(left-hand) electric polarization and positive (negative) helicity. It also finds a > 0 for the 
Alfven mode in /? = 1 plasmas. Moreover, when the wave tends to more oblique propagation, 
the ion cross-helicity of the slow magnetosonic wave aci ^ — 1 . 

Fig. (4) shows that the Alfven and slow modes are circularly polarized {\6Ex\ — \SEy\) 
waves at cj < Ud in the high-/? plasmas, where the Alfven (slow) mode exhibits the right- 
hand (left-hand) electric polarization and positive (negative) helicity. These two modes also 
have the same ion cross-helicity distribution. Note that three modes have no interaction 
point at the high oblique propagation as shown in Panel (cl). 

It needs to note that the electric held polarizations also strongly depend on the ratio 
of the electron to ion temperature Te/Tj (Eq. fl2U]) h The electric held polarizations with 
diherent Tg/Tj are presented in Fig. (5), where 6 = 45° and Te/Tj = 0, 1, and cxo. It 
shows that the parallel polarization \6Ez/6Ey\ decreases obviously with decreasing Tg/Tj. 
The transverse polarization \6Ex/6Ey\ is slightly ahected by Tg/Tj for the long-wavelength 
Agfc <C 1 waves, but not for the long-wavelength fast mode in the high-/? plasmas or for the 
long-wavelength slow mode in the low-/? plasmas. To understand qualitatively above results, 
the complete expressions Eq. (|20D can reduce to 


SE^ ^ + Tg (o;^ - V^k^) + Vjk^k^,] 

i6Ey (cj^ ~ ^A^z) 

5E, ^ + 

i6Ey k± oJOJd 


(29) 


where the long-wavelength (Ag/c 1) and very low-frequency (Ai 2 — 1) conditions are used, 
and the smaller terms of the order of Q are neglected. Since Tj = 1 / (1 -|- ygEg/q^Tj) and Tg = 

— Tjj, \6Ez/ {i6Ey)\ decreases with decreasing Tg/Tj. When 00 ~ VAkz, SE^/ {iSEy) ~ 

—Vj^k^Xikz/ (ca^ — V^kl) is independent on Tg/Tj. We can also hnd dE^/ {iSEy) ~ —T^u/ud 
corresponding to the fast wave up ~ Vpk in /? 3> 1 plasmas and SE^/ {iSEy) ~ T^V^k'^/uud 
corresponding to the slow wave ojs ~ Vpkz in /? <^ 1 plasmas, which indicate \5Ex/ {i5Ey)\ 
decreasing with smaller Tg/Tj. Besides, in Fig. (5) the electric polarizations \5Ex/5Ey\ and 
\5Ez/5Ey \ in Tg/Tj 7 ^ 0 plasmas increase continuously as the slow magnetosonic wave changes 
to the ion cyclotron wave, while both polarizations are nearly unchanged in Tg/Tj = 0 
plasmas. Note that the main characters in Fig. (5) still appear in the electric polarization 
distributions of the oblique waves with 9 7 ^ 45°. 

When the phase relation of the electric polarizations changes, a peak or a valley can occur 
in \5Ez:/5Ey\ and \5Ez/5Ey\ distributions. For the fast mode in Fig. (5), the phase relation 
of 5Ez/5Ey changes in T^/Ti = 0 plasmas, while it is unchanged in Tg/Tj 7 ^ 0 plasmas. For 
the Alfven mode, two transition points arise in the phase relation of SE^/SEy in T^/Ti 7 ^ 0 
plasmas; however, in the cold electron {T^/Ti = 0 ) plasmas, the transition at the smaller Ajfc 
disappears in /? <C 1 plasmas, or two transitions are both missing in /? > 1 plasmas. 

For near-perpendicular Alfven wave. Fig. ( 6 ) shows that the parallel polarization 
\6Ez/SEy\ increases (decreases) with Tg/Tj as Aj/c < 10^ (Aj/c > 10^). At 9 = 89.99°, the 
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FIG. 5: The sensitivity of the electric polarizations of the three modes on the temperature ratio 
Tg/Ti, where the waves propagate at 9 = 45° in the plasmas with different /3 environments: low-/3 
(/3 = 10“^), /3 = 1, and high-/3 (/3 = 10^). The solid, dashed and dotted lines represent the cold 
electron T^/Ti = 0, equal ion and electron temperature Tg/Tj = 1, and cold ion Tf./Ti = oo cases, 
respectively. 


transverse polarization \6Ex/5Ey \ decreases with increasing Te/Tj for the kinetic-scale Alfven 
waves (Ajfc > 1). These resnlts indicate the important role of the electron temperatnre Te on 
the kinetic-scale Alfven waves BUS- Moreover, there is no transition of the phase relation 
of 6Ex/5Ey at 6 = 89.99° case. The reason is that the wave freqnency u is smaller than 
the ion cyclotron freqnency ujd at 9 = 89.99°, which cannot satisfy the freqnency condition 
oj > ujci for the changing of the phase relation of transverse electric polarization (lll |. 
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FIG. 6: The dependence of the electric polarizations of high oblique Alfven waves on the temper¬ 
ature ratio Te/Ti, where two angles 6 = 85° and 89.99° are considered. The solid, dashed and 
dotted lines represent the cold electron Tg/Tj = 0, equal ion and electron temperature Tg/Tj = 1, 
and cold ion Tg/Tj = oo cases, respectively. 


IV. SUMMARY 

In this study ions and electrons are treated separately in comparison with one fluid 
element {U = ^ rriaVa/ '^a) method adopted in previous studies n This method is 
helpful to obtain the linear eigenfunctions including the ion and electron velocities as well 
as the ion and electron cross-helicities. 

It found that the fast and Alfven modes are nearly linearly polarized at the very low- 
frequency uj <C 10~^ojci, and circularly polarized at 10~‘^ujci < oj < ujd at the near-parallel 
propagation in the low-/3 plasmas. Two circularly polarized modes become the fast and slow 
modes in a narrow frequency regime O.lcjci < oo < ood in /? = 1 plasmas; they are the Alfven 
and slow modes m. (3 = 10^ plasmas. To the ion cross-helicity aci of the long-wavelength 
slow mode, aci — 0 in the low-/? plasmas, aci —)■ — 1 as 0 —)■ 90° in /? = 1 plasmas, and 
aci — —1 in the high-/? plasmas. It also found that the negative magnetic-helicity a of the 
Alfven mode can occur at the small or moderate angles in the low-/? plasmas, while a > 0 
arises always at the high oblique angle in the low-/? plasmas or at the general angle in /? > 1 
plasmas. 

Our results exhibited the sensitivity of the electric polarizations on the temperature 
ratio Te/Ti. The parallel polarization \5Ez/5Ey\ decreases with Te/Ti as \ik < 1. The 
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transverse polarizations \6Ex/5Ey\ also decreases with Tg/Tj for the long-wavelength fast 
mode in the high-/9 plasmas, or for the long-wavelength slow mode in the low-/? plasmas, 
while \6Ex/5Ey\ at other long-wavelength cases are slightly affected by Te/Tj. Fnrthermore, 
the phase relation of 5Ex/5Ey of the Alfven mode will change in Te/Tj ^ 0 plasmas, bnt this 
change can disappear in the cold electron T^/Ti = 0 plasmas. For the fast mode, the phase 
relation of 6Ez/5Ey changes in Tg/Tj = 0 plasmas, while the nnchanged phase relation arises 
in Te/Ti 7 ^ 0 plasmas. 

We have also presented the approximate dispersion relations in the near-perpendicnlar 
propagation, low-/?, and high-/? limits. These approximations can describe nicely the exact 
dispersion relations of the three modes given by Eq. ffTbjl . It notes that the condition of 
V/l/c^ -C 1 nsed in the stndy leads to the neglecting of the displacement cnrrent. However, 
this assnmption is broken near the wave cntoff position which resnlts in the validity condition 
of a; <C Uci (1 + (? IVa) 0 . Also, the displacement cnrrent may be important in prodncing 


the parallel electric held of the low-freqnency Alfven mode [1^. Therefore, a comprehensive 


stndy inclnding the effect of the displacement cnrrent is needed. 

Lastly, two-hnid model neglects the kinetic wave-particle interaction effects, snch as Lan- 
dan damping and ion (electron) cyclotron resonance damping, which can only be captnred 
by the kinetic model. These kinetic effects can strongly affect the wave dispersion relation 
and polarization properties. For example, the wave dispersion relation of the kinetic Alfven 
wave is depressed at ion scales in the high-/? plasmas where there can be the heavy Landan 
damping The two models also result in different phase relation between two electric 
components. However, since it is difficult to identify clearly all modes from the full kinetic 
theory, the two-hnid theory can be a useful guide to discard the modes in the kinetic theory. 
Our complete expressions can be conveniently used to compare with the results of the kinetic 
model. 


Appendix A: Approximate dispersion relations in different limits 
1. Near-perpendicnlar propagation limit 

At the near-perpendicular propagation limit k‘l/k‘]_ <C 1, the cubic equation for 
can reduce to a quadratic equation 

Aobiii^^ — TobiiW^ -|- Oobii = 0, (Al) 

where 

Aobii = 1 + Xlk'^ + X^ikl + (1 + Xlk'^Y 13 + kllk\ 

Tobu = {l + 2f3 + p^k‘^)Vlkl 

C-obii = l3V^kt, 

which contains the dispersion relation of the Alfven mode {ua = i^+) and slow mode 
(cjs = a;_). 
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1 + 213 + 


4 = vikt 


2 (1 + + Xfkl + (1 + XlkY P + kl/k‘^) 


l±Jl-A(3 


1 + Xlk-^ + Xfkl + (1 + Xlk^f (3 + kl/k^ 


(l + 2/? + p2fc2)2 

On the other hand, the fast mode can be obtained by hrst two terms in Eq. ffTHll . 


(A2) 


Up = VAk 


^ 1 + A2F + Xfkl + (1 + Xlk^f 13 + kl/k^ 


1 + Xlk^ 


(A3) 


2. Low-/3 (/3 ^ 1) limit 


In low-/? plasmas with /9 1, the freqnency of the slow mode ojs ~ Vpk^ is mnch smaller 

than that of the fast mode up ~ and Alfven mode ua ~ So that the last two 

terms in Eq. flT^ control the dispersion relation of the slow mode 


Vxk, 

^yl + p‘^k3 


(A4) 


On the other hand, the fast and Alfven modes are controlled by following qnadratic eqnation 


~ B\ow^‘^ + C'low — 0, (As) 

where 

Alow = (l + ) 

Slow = + ^ 2 /^^) ^4^^) 

Olow = {l+P^e)Vjekl 

which yields the dispersion relations for the fast mode ( 0 ;^ = a;+) and Alfven mode 
{ua = UJ-) 


u, = 


Vlk^ 


1 + xie + (1 + xik^f /j + Xfki + kijk^ 

2{1+ Xlk^f 


1 ± 


A{l + Xlk‘^f{l + p^k^) 


p2 


\ (1 + A2fc2 + (1 + xik^fp + xfki + ki/k^y 


k^ 


(A6) 


3. High-/3 (/3 ^ 1) limit 

In high-/? plasmas with /3 ^ 1, the freqnency of the fast mode up ^ Vpk is mnch larger 
than that of the Alfven and slow modes oja,s ~ The fast mode are dominant by 

the hrst two terms in Eq. fflSj) . whereas the Alfven and slow modes are dominant by the 
qnadratic eqnation shown in Eq. flAip . Therefore, the wave dispersion relations of the three 
modes are the same as that given by Eqs. flA2p and flASp . 
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